Magnetoelectric effects in heavy-fermion superconductors without inversion symmetry 
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We investigate effects of strong electron correlation on magnetoelectric transport phenomena in 
noncentrosymmetric superconductors with particular emphasis on its application to the recently 
discovered heavy-fermion superconductor CePtaSi. Taking into account electron correlation effects 
in a formally exact way, we obtain the expression of the magnetoelectric coefficient for the Zeeman- 
field-induced paramagnetic supercurrent, of which the existence was predicted more than a decade 
ago. It is found that in contrast to the usual Meissner current, which is much reduced by the mass 
renormalization factor in the heavy-fermion state, the paramagnetic supercurrent is not affected by 
the Fermi liquid effect. This result implies that the experimental observation of the magnetoelectric 
effect is more feasible in heavy-fermion systems than that in conventional metals with moderate 
effective mass. 
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I. INTRODUCTION 

It has been discussed for decades that metallic systems 
with noncentrosymmetric crystal structure may exhibit 
non-trivial magnetoelectric effects The existence 

of an asymmetric potential gradient VU due to the non- 
centrosymmetric structure gives rise to the spin-orbit in- 
teraction (p x VV) • ex, which breaks the individual in- 
version and spin rotation symmetry. As a result, the 
charge or energy current operator may couple to the spin 
density operator. In this context, current flows induced 
by an applied magnetic field, and current-flow-driven 
magnetization have been investigated extensively both in 
normal metalsi^ and superconductors^^ In particular, 
Edelstein predicted the remarkable magnetoelectric ef- 
fect in superconducting states; i.e. in noncentrosymmet- 
ric superconductors, the Zeeman field induces a super- 
current, and vice versa, the supercurrent flow induces a 
magnetization,^ Later, the former effect is elegantly re- 
formulated by Yip in terms of the "van Vleck" contribu- 
tion which stems from the inversion-symmetry-breaking 
spin-orbit interaction. 5 Since a static magnetic field can 
not induce dissipative current flows, the Zeeman-energy- 
induced current should vanish in the normal state. How- 
ever, in the superconducting state, the existence of the 
paramagnetic supercurrent is not forbidden in the ab- 
sence of the inversion symmetry. The recent discovery 
of superconducting materials without inversion symme- 
try such as CePtaSi, Ulr, and Cd2Re2C>7 stimulates the 
renewed interest in this i ssuei £iL2i2ii2iiiii£ Under an ap- 
plied magnetic field, the Meissner diamagnetic supercur- 
rent in addition to the Zeeman-field-induced paramag- 
netic supercurrent should exist. Thus, it is important for 
the experimental observation of this effect to discriminate 
between these two supercurrents. 

In the present paper, we would like to investigate the 
Fermi liquid corrections to this Edelstein's magnetoelec- 
tric effect, which may be important for the application 
to heavy-fermion superconductors such as CePtaSi and 
Ulr. We obtain the formula for the magnetoelectric effect 
coefficient taking into account Fermi liquid corrections 
exactly. The most important finding is that the Zeeman- 



energy-induced paramagnetic supercurrent is not at all 
affected by electron correlation effects provided that fer- 
romagnetic spin fluctuation is not developed, in contrast 
with the diamagnetic Meissner current of which the mag- 
nitude is much reduced by the mass renormalization ef- 
fect. This result implies that the experimental detection 
of the paramagnetic supercurrent in heavy-fermion su- 
perconductors may be more feasible than that in weakly 
correlated metals. 

The organization of this paper is as follows. In Sec. II, 
we present the basic formulation of the Fermi liquid the- 
ory for a model system without inversion symmetry. We 
would like to make a brief comment on the superconduct- 
ing state realized in CePtaSi in Sec. III. In Sec. IV, the ex- 
act formula of the magnetoelectric coefficient is obtained. 
In Sec.V, the implication for the experimental observa- 
tion of this effect is discussed. Summary and discussion 
are given in the last section. 



II. MODEL AND ANALYSIS BASED ON THE 
FERMI LIQUID THEORY 

As a simplest model which realizes the broken inver- 
sion symmetry, we consider an interacting electron sys- 
tem with the Rashba spin-orbit interaction^ 3 The Hamil- 
tonian is given by, 

+ u ^2 n H n li> (!) 

i 

where $ (c ap ) is the creation (annihilation) operator for 

The number 
The second 



an electron with spin a and momentum p. 
density operator at the site i, n a i 
term of Q is the Rashba spin-orbit interaction which 
incorporates the broken inversion symmetry. Here, the 
unit vector parallel to the asymmetric potential gradient 
is given by ft = (0, 0, 1). This system is considered to be 
a model of CePtaSi, with which we are mainly concerned 
in this paper. The /-electron of CePtaSi is in the Tj 
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Kramers doublet stated Expanding the doublet in 
terms of the s z = 1/2 and —1/2 basis, we found that 
the Rashba spin-orbit interaction term expressed in term 
of the I\ basis has the matrix structure given above up 
to a constant factor which can be absorbed into the re- 
definition of the coupling constant a. Thus, in the case 
of CePtaSi, the spin index a in represents the 
Kramers doublet. 

In the following, we do not specify the pairing mecha- 
nism of superconductivity, but assume that the supercon- 
ducting state is realized by an effective pairing interaction 
with an angular momentum I > 1 which may stem from 
the on-site Coulomb interaction in or may have any 
other origin not included in the Hamiltonian (JJJ. Then, 
we can analyze electron correlation effects on this super- 
conducting state in a formally exact way by using the 
superconducting Fermi liquid theory^*i^& 

In the conventional Nambu representation^ the in- 
verse of the single-particle Green's function is defined as, 



a-u n \ = ( i£ ™- h (p) -Mp) 

{P1 \ -At(p) ie n + H t (-p) 
where p = (p, is), and, 

H( P ) = H (p) + £(p), 

Ho = e p - (j, + a(p x ft) ■ a - p, B cr x H x 



(2) 



(3) 



(4) 



with p chemical potential. Here, to discuss the magneto- 
electric effect, we take into account the Zeeman magnetic 
field in the x direction, H x . For simplicity, we assume 
that the (/-value is equal to 2. The self-energy matrix £ 
consists of both diagonal and off-diagonal components, 



( £ TT (» £ u (» 
- - 1 E iT (p) E u (p) 



(5) 



We can easily see from the symmetry argument that 
under the applied in-plane magnetic field, Xff(p) = 

s ub) = and s it(^* £ ) = ^n(p,- i£ )- Tnc 

superconducting gap function is {A(p)} a p = A a p(p) 
(a, 0=1,1). 

ie n — H(p) and «e n +i?*(— p) in G~ x (p) are diagonalized 
by the unitary transformation A{p)G~ 1 {p)A^ (p) with, 



A[P) ~ ^ U^(-p) 



U(p) 



=-( 1 



it- 
1 



(6) 



(7) 



where t± = t x ± it y , and t x , i y are, respectively, the x 
and y components of the unit vector t p = t p /\t p \ with 

t p (ie) = ( Px + hmZ u ,p y + ^Re£ n - ^^,0). (8) 



As seen from eqs.QjQ, and JSJ), the main effect of 
the off-diagonal self-energy is to renormalize the 
Rashba interaction term, replacing the momentum p in 
the Rashba term with the vector t p . Since the on-site 
Coulomb interaction does not change the symmetry of 
the system, the off-diagonal self-energy should satisfy the 
following condition in the absence of magnetic fields; 



ReS TJ ,(p x , -p y ) = -BBE n (p x ,p y ), 
lraE n (-p X) p y ) = -ImE n (p x ,Py), 



(9) 



(10) 



and ReEjj (ImS-^) is an even function of p x (p y ). 

In the normal state, the single-particle excitation en- 
ergy e* T for the quasi-particle with the helicity t = ±1 

is given by the solution of the equation z — H{p, z) = 0, 
which is, in the diagonalized representation, 

e* pr + M - s p - ra\t p (s; T )\ - ReE® s* pT ) = 0. (11) 

The gap functions A(p) and A^(p) in G^ip) 
are also diagonalized by the unitary transformation 
A{p)G^ 1 {p)A^ (p) provided that the gap function has the 
following structure, 



A(p) = A s (p)ia y + A t (p)(i p x n) 



(12) 



Here A s (p) and A t (p) are even functions of momentum p. 
This means that the spin singlet and triplet component 
is mixed in the diagonalized basis labeled by r = ±1, and 



the d vector of the triplet component is L 



i. 3 ' 18 In the 



case that A(p) is not expressed as ea. <jT2T) . H (p) and A(p) 
can not be diagonalized simultaneously, and the non-zero 
off-diagonal components of A(p) which correspond to the 
Cooper pairing between the different Fermi surfaces in- 
duce pair-breaking effects, resulting in the decrease of 
the transition temperature T c . Thus, the highest transi- 
tion temperature is achieved by the gap function given 
by ea. l|12ll ii2i The realization of the gap function (|12fl in 
the case with no inversion center is also elucidated by the 
group theoretical argument^ 

Taking the inverse of eq.©, we have, 



5(p) 



( G{p) F(p) 
\P(p) -G\-p) 



where 



G(p) 



E 

r=±l 



1 + T(i p x n) 



-Grip), 



F iP) = 1^ 9 i<?yFrip), 



T = ±l 



and, 



Grip) 



z P rjie + e* pT ) 
{is + ijsgne) 2 - E 2 ' 



(13) 



(14) 



(15) 



(16) 
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F T (P) 



E 2 ' 



(is + ijsgne) 2 
Here the mass renormalization factor is, 
9ReE(p) 



(17) 



1 - 



d{ie) 



aimS u (p) f 9RcS Ti (p) 
' ' * 8(ie) + y d(ie) 



ie=E pr 



and 7 is the quasi-particle damping. The single-particle 
excitation energy is E pT = ^Je* 2 + A 2 , (p) with A pr = 
(A.(p)+rA t (p)). 

The superconducting gap function A Q( g and the tran- 
sition temperature are determined by the self-consistent 
gap equation, 



(19) 



n.p 



where we have introduced the four-point vertex func- 
tion matrix {f"^ (PiP')}fS which is diagrammatically ex- 
pressed as shown in FIG.l. We expand the four-point 
vertex in the particle-particle channel {T a P (p,p')}js m 
terms of the basis of the irreducible representations of 
the point group, and consider a component {T"^ (p,p')}jS 
which corresponds to the pairing state giving highest T c . 
This pairing interaction consists of the spin singlet and 
triplet channel: 

{^Z P {p,p')}iS = r 's{j>,p')i{oy) a f3i{a v ) 1 8 
+T t (p,p')(i p x n) ■ {aio- y ) a p(tp> x n) • (Si(j y ) lS {20) 

The symmetries of T s (p,p / ) and T t (p,p') in the momen- 
tum space are characterized by the same irreducible rep- 
resentation, and as a result, A s (p) and A t (p) in ca. ltT^|) 
have the same symmetry in the momentum space. The 
realized superconducting state is the mixture of the spin 
singlet and triplet statesAi^ In this case, the possible 
pairing state is s +p or d + / or g + h, and so forth. 



III. A COMMENT ON THE 
SUPERCONDUCTING STATE REALIZED IN 
CePt 3 Si : POSSIBILITY OF AN UNUSUAL 
COHERENCE EFFECT 

Here, we would like to make a brief remark about the 
pairing state realized in CePtaSi. The NMR measure- 
ment carried out by Yogi et al. shows the existence of 
the coherence peak of \fT\T just below T c , indicating 
the full-gap state without nodes^ On the other hand, 
the recent experiment on the thermal transport done 
by Izawa et al. supports the existence of line nodes of 
the superconducting gap^ A possible resolution of this 
discrepancy is that the line nodes of the superconduct- 
ing gap are generated accidentally at the magnetic Zone 



boundary which emerges as a result of the antiferromag- 
nctic phase transition at Tn = 2.2 K, and crosses the 
Fermi surface. For such accidental nodes without the 
sign change of the superconducting gap function, the co- 
herence factor of \jT\T does not vanish, resulting in the 
enhancement of the coherence peak just below T c . In this 
case, a plausible candidate for the pairing state is the s+p 
wave state. An important point which we would like to 
stress here is that even when the superconducting state is 
dominated by the p wave pairing; i.e. A s (p) <C A t (p), the 
coherence factor which enters into \jT\T does not van- 
ish. This contrasts with the case of the usual p wave state 
realized in centrosymmetric superconductors, where the 
coherence factor of \jT\T disappears. This is understood 
as follows. For simplicity, we ignore electron correlation 
effects. Then, in noncentrosymmetric superconductors, 
the nuclear relaxation rate is given by^£ 

km llm[TVV{Tr[^G(p, e m + u n )?-G(p', e m )\ 

e m p,p' 

-Tr[— F(p,e m + w„)y F(p , £ m )]}\iu n ^+iS 



de 



1 



2tt 2Tcosli ^ 



{[iV n ( £ )] 2 + [7V a (s) 2 ]}, 



(21) 



with N n (e) and N a (e) defined by the retarded Green's 
functions as, 



JV„(£),-^^ImG r K ( P , £ ), 



N a (e) 



p r — ± 



(22) 



(23) 



p r — ± 



The expression of 1 /T\T (|21|l does not rely on the phase 

factor tpXn of the triplet component of the gap function 
(|12|l . The second term of the right-hand side of (|21|l gives 
the non-zero contribution from the coherence factor, as 
in the case of conventional s wave superconductors. This 
property enhances the coherence peak of 1/T\T promi- 
nently. It may be important to take into account this 
unusual coherence effect for clarification of the origin of 
the notable coherence peak of \jT\T in CePtaSL 12 It is 
also intriguing to explore the unusual coherence effect on 
other response functions, such as the ultrasonic attenua- 
tion. We would like to address this issue elsewhere. 



IV. MANY-BODY EFFECTS ON THE 
MAGNETOELECTRIC TRANSPORT IN THE 
SUPERCONDUCTING STATE 

In this section, we consider electron correlation effects 
on the magnetoelectric transport in the superconduct- 
ing state first found by Edelstein and later discussed in 
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-P' ,5 



p ,a 



-P,P 



FIG. 1: Four point vertex in the particle-particle channel 



where the three-point vertex function is 
A£_(p)=* B (l 



1 aReSti, L 9ImE T , 

+ — Q ' ■ (32) 



The last two terms of (I31[) . which emerge as a result 
of Fermi liquid corrections, can be rewritten by using 
the Luttinger- Ward's identity generalized to the super- 
conducting state. The Luttinger- Ward's identity in the 
normal state readsri* 2 ^ 



detail by Yip; i.e. the emergence of the paramagnetic su- 
percurrent induced by the Zeeman magnetic field in the 
direction normal to the n vector£& 

We consider the charge current flowing in the y direc- 
tion which is defined as, 



with 



Vr 



Oy 



V0y(p) 

-^(-P) 



(24) 



(25) 



and Vq v = d Py £ p + a(n x a). 

The transport coefficient which characterizes the Edel- 
stein's magnetoelectric effect is given by, 



_ 8Jy 



dH x 

1 

2 



H x =0 



n,p 



dH x 



Following ref.3, to simplify the expression of JC yx we 
use the Ward's identity for the current vertex, 

'' (:ip] - G(p)Vy(^G(p) + F(p)V t y (-p)F\p,- i e) 



d Py 



-R(p), 



where V y {p) — i>o y + 9£(p) / 'dp y , and 

rip) 



r{p) 



R{p) = 



2a\t r 



r*ip,-ie) 
- G+G- + F+K 



(27) 



(28) 



Aff_(p) =*«(<* + Q „ U )-f 



dPy 8py 

Then, from lfP3)l. and (J27J), we obtain, 



i + A c + y _(p), (29) 
(30) 



K 



-2T^Tr 



Fip)vl y {-p)F\p) 



2n B aTj2F+Fj x A s *_ip) 

7 

T^Tr 



n.p 



' dt 


dG ' 


-T^Tr 

n,p 


' dt 


dG 


DPy 


dH x 


dH x 


d Py_ 



(31) 



(33) 



This relation is obtained by differentiating all closed 
linked diagrams with respect to p M . 21 ' 22 In the super- 
conducting state, a similar analysis leads, 

n,p ** n? p ^ 

+ rVTr[A^]=0. (34) 
dp,, 

Differentiating eq.J2U) with respect to H x , and integrat- 
ing over p y by parts, we have, 



T^Tr 



' dt 


dG ' 


-T^Tr 

n,p 


' dt 


dG~ 


d P y 


8H X 


dH x 


d Py_ 



-T^Tr 



n,p 



-T^Tr 



aAt of 

dp y dH x 

dA <9#t 
dp y dH x 



n , p 



+ T^Tr 



<9At dF 
dH x dp y 

dA dP^ 
dH x dp y 



(35) 



We see from eas. (l31|l and (|35|l that JC yx vanishes exactly 
in the normal state, as is consistent with the thermo- 
dynamic argument that a static magnetic field can not 
induce non-equilibrium current flowsi^ The right-hand 
side of eq. (|35|l consists of the terms which have the form 
J2 n p G T F T r . The ratio of the contributions from these 
terms to those from other terms of ea. (|3I|l is of order 
A/Ep, and thus we can neglect the last two terms of l|31|) 
approximately. Then, the magnetoelectric coefficient is 
expressed as, 



fcyx 



z 2 A 2 

pr pr 



p r— ±1 



pT 



ch 



-2 E pT 
IT 



th- 



2T 



„ z p+ z p ^ A p+ A p _ 
+2a ^ E*-El_ 



p 



2T 



th%t 



th 



2T 



E. 



p+ 



E,, 



where 



dH x [is dH. 



1 1 



K x (p) 

i x A%%m 



MB C-Hx 



5 



It is noted that the vertex corrections due to electron 
correlation, A sx (p), which appear in the above expression 
(|36|) is nothing but the vertex corrections to the uniform 
spin susceptibility, 



Xxx 



,2 \ \ ^ 



P T = ±l P T 



*2 ch -2E£Z. 

2T 



A* th- 



2T 



E pT 



iyK X (p) 



-4E 



f* th p+ 
P+ 2T 



2E„ 



-* th Ep ~ 

2_E„_ 



a\t 7 



:K X -(P)- (38) 



Ea. (|38[l is easily obtained by differentiating the 
x component of the total magnetization S x — 
HbTJ^u P Tr[ca;G(p)] with respect to H x . Note that the 
first term of the right-hand side of i|38|) is the Pauli para- 
magnetic contribution and the second one is the "van 
Vleck" term which arises from excitations between spin- 
orbit split two bands. Generally, in heavy- fermion sys- 
tems, the magnitude of the uniform spin susceptibility 
is enhanced by the vertex corrections A sx (p). In typical 
heavy-fermion systems including CePtaSi, the Wilson ra- 
tio R w = Tx/C/(Txo/C ) ~ 2, which implies that the 
vertex corrections A sx is approximately of order the mass 
enhancement factor l/z pT & Therefore, in ea. l|36[) effects 
of the vertex corrections and the mass renormalization 
factors z pT cancel with each other. This cancellation 
holds as long as there is no strong ferromagnetic spin 
fluctuation which increases notably the magnitudes of 
the vertex corrections A sx . Another important feature of 
eq. (|36[) is the absence of the backflow term of the charge 
current, which usually exists in the non-equilibrium cur- 
rent flow. (See the discussion on the usual Meissner cur- 
rent in the next section.) This is related to the fact that 
the current induced by a static magnetic field is a dissi- 
pationless equilibrium flow. As a result, the Fermi liquid 
corrections do not exist in this magnetoelectric coefficient 
for heavy fermion superconductors, provided that there is 
no ferromagnetic fluctuation. This is one of the main re- 
sult of this paper. In terms of the Kubo formula, the ab- 
sence of electron correlation effects for JC yx is understood 
as follows. K-y X is given by the correlation function of the 
current and spin density operators. The spin density ver- 
tex is renormalized by electron correlation in the opposite 
way to the current vertex, resulting in the cancellation 
of the mass renormalization factors. The important im- 
plication of this result is that the Zeeman-field-induced 
paramagnetic supercurrent is not suppressed by electron 
correlation effects in contrast with the usual diamagnetic 
supercurrent of which the magnitude is much reduced by 
the large mass enhancement in heavy fermion systems. 
This property may make the experimental observation of 
the magnetoelectric effect easier, as discussed in the next 
section. 



V. IMPLICATIONS FOR EXPERIMENTAL 
OBSERVATIONS 

On the basis of the formula (|36|) , we would like to dis- 
cuss how the Zeeman-field-induced paramagnetic super- 
current is experimentally observed. 

When the magnetic field is applied in the x direction, 
the London equation is modified to, 



Js 



47rA ; 



:A + K,y X (n x H x ), 



(39) 



where A is the vector potential. Since the applied mag- 
netic field always induces both the diamagnetic and 
the paramagnetic supercurrent in the system with the 
Rashba spin-orbit interaction, it is important for the ex- 
perimental observation of this effect to discriminate be- 
tween these two supercurrents. If one measures currents 
simply attaching leads to the sample and applying an 
in-plane magnetic field, the Zeeman-field-induced para- 
magnetic supercurrent can not flow in the sample because 
it generates the Joule heat in the normal metal leads' 
In this case, the magnetoelectric effect is canceled with 
the non-zero phase gradient of the superconducting order 
parameter in the equilibrium state. 

Here, to highlight the observation of the paramagnetic 
supercurrent, we consider an experimental setup com- 
posed of two superconducting samples joined together as 
depicted in FIG. 2. The n vectors of the sample A and B 
are, respectively, given by fix = (0, 0, 1), and % = — ua- 
The joined surface at the junction is normal to the n 
vectors. The applied magnetic field in the x direction 
H = (H, 0, 0) gives rise to the paramagnetic supercurrent 
in the sample A (B) in the direction ua X H (-% x H). 
Then, electrons accumulated in the right (left) edge of the 
sample A (B) transfer to the right (left) edge of the sam- 
ple B (A) to decrease the chemical potential difference 
between the samples A and B, resulting in supercurrent 
flows circulating in the system. The applied magnetic 
field also induces the Meissner diamagnetic supercurrent. 
As explained below, the paramagnetic contribution can 
be discriminated from the diamagnetic current by using 
the Volovik effect^ 

Let us first consider the coefficient of the diamagnetic 
Meissner supercurrent with the Fermi liquid corrections 
at zero temperature, which is equal to the Drude weight 
in the normal stated 



4ttA 2 



v£J£5(p-e*), 



(40) 



where the quasi-particle velocity is v^* = de^/dp^, and 
the charge current is 



(41) 



Here f pT ,p'T' is the interaction between two quasiparti- 
cles. The second term of (I41|l is the backflow term. In 
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heavy fermion systems, the mass renormalization factor 
z pT and the backflow term in the current J* T give rise 
to the pronounced reduction of the Meissner coefficient. 
For example, in CePtaSi, z pT is estimated as ~ l/lOOp^iS 
If we assume the spherical Fermi surface, ea. (l40|) reduces 
to (e 2 /c)v* F n a /pF- Here v* F is the renormalized Fermi 
velocity, and n s the superfluid density, p F the Fermi mo- 
mentum. It should be notified that n s is renormalized 
by the backflow effect of J* T , and is not equal to the car- 
rier density even at zero temperature. In particular, for 
heavy fermion systems in which Umklapp scattering is 
expected to be strong, n s is smaller than the carrier den- 
sity. In contrast to the diamagnetic Meissner current, 
the Zeeman-field-induced paramagnetic supercurrent is 
not influenced by the many-body effects described above, 
as discussed in the previous section. This difference of 
electron correlation effects between the diamagnetic and 
paramagnetic supercurrents can be utilized to detect the 
magnetoelectric effect. 



+ 



n A 



B 

•* — y 



i 



n B 



FIG. 2: An experimental setup for the detection of the 
Zeeman-field-induced supercurrent. The n vectors of the two 
superconducting samples A and B (depicted by the gray ar- 
rows) are aligned in the directions (0, 0, 1) and (0, 0,-1), re- 
spectively. The in-plane magnetic field H is applied in the 
x direction. The paramagnetic supercurrent circulates in the 
system as depicted by the black thin arrows. 



If the superconducting gap has a nodal structure as 
is often realized in some heavy-fermion superconductors, 
the existence of the paramagnetic supercurrent is indi- 
rectly observed through the Volovik effect on the single 
particle density of states under the applied in-plane mag- 
netic field for H c i < H < H C 2. Here H c \ and H C 2 are, re- 
spectively, the lower and upper critical field. In fact, the 
recent thermal transport measurements for CePt 3 Si sup- 
ports the existence of the line node in the superconduct- 
ing state of this system. 23 Applying the semi-classical ap- 
proximation based upon the Doppler shift effect^ and 
assuming a spherical Fermi surface, we calculate the lo- 
cal density of states from the modified London equation 
(|39|l . In the calculation of JC yx , we use the fact that the 

vertex corrections A*^ ^ is appropriately approxiamted 

as 

typical heavy-ferimon systems as discussed 
in Sec. IV., and expand ea. (|36fl in terms of app/Ep up 



to the lowest order. The result at T = is 



•JAoc(O) 



H, 



<::2 



± 



H C 2 TT^EfUs 



• (42) 



Here $o = hc/(2e), and no is the density of electrons. 
Ep is the unrenormalized Fermi energy. The first term 
of the right-hand side of l|42[l is due to the usual Volovik 
effect, and the second term linearly proportional to the 
applied magnetic field stems from the Zeeman- energy- 
induced paramagnetic supercurrent. Thus, the magnetic- 
field-dependence distinguishes between the paramagnetic 
and diamagnetic currents. The above behavior of the lo- 
cal density of states may be observed by the measurement 
of the specific heat coefficient or the thermal conductivity 
in sufficiently low field regions i2&SL22i2S In the above ex- 
pression of 8Di oc (0), it is seen that the conventional dia- 
magnetic contribution is suppressed by the mass renor- 
malization factor z pT which appears through v* F , whereas 
the magnetoelectric contribution is not affected by this 
correlation effect. It is also noted that the carrier density 
which enters into the paramagnetic term is not n s but 
equal to the electron density no- This is due to the ab- 
sence of the backflow term in the Zeeman-energy-induced 
supercurrent. As mentioned before, n s is affected by the 
backflow term. For simplicity, we assume that n s w uq 
for a while. 

In the case of CePtaSi, according to the measurement 
of H C 2, the coherence length £ ~ 8.1 x 10~ 7 chip It is a 
bit difficult to estimate the renormalized Fermi velocity 
from experimental measurements. Bauer et al. obtained 
v* F ~ 5.29 x 10 5 cm/s from the data of dH&ldT and 
the specific heat coefficient, assuming a spherical Fermi 
surfaced This value of v* F is almost of the same order as 
that obtained by combining the unrenormalized Fermi 
velocity computed from the LDA method and the mass 
enhancement factor z _1 ~ 100 estimated from the spe- 
cific heat measurement Ai£ According to the LDA band 
calculationSfiS* 3 ^ the spin-orbit splitting is not so small 
compared to the Fermi energy, and may be approximated 
as app/Ep ~ 0.1. Then, for CePtaSi, we have, 



*Aoc(o) 



1.0 x 10- 24 ,/-^- ± 0.48 x 10 - 24 -^ 

V Hc2 Hc2 



.(43) 



It is remarkable that the contribution from the param- 
agnetic supercurrent (the second term of l|43l) ) is compa- 
rable to that from the Meissner supercurrent (the first 
term of l@3Jl). It should be stressed that the feasibility 
of the experimental observation of the Edelstein's mag- 
netoelectric effect is due to the large mass enhancement 
in the heavy fermion system, which suppresses strongly 
the Meissner supercurrent, but in contrast, does not 
affect the Zeeman-field-induced paramagnetic supercur- 
rent. Moreover, if the superfluid density n s is reduced 
by the backflow effect, the Meissner term of (|42|) is more 
suppressed compared with the paramagnetic term, and 
thus the observation of the magnetoelectric effect may 
become easier . 
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VI. SUMMARY AND DISCUSSION 

We have investigated electron correlation effects on 
the magnetoelectric transport phenomena in supercon- 
ductors without inversion symmetry, ft is found that, 
in contrast to the Meissner diamagnetic supercurrent 
which is much reduced by the mass enhancement factor in 
the absence of translational symmetry, the Zceman-field- 
induced paramagnetic supercurrent is not affected by the 
strong electron correlation provided that ferromagnetic 
fluctuation is not developed. Because of this remarkable 
property, the experimental detection of the magnetoelec- 
tric effect may be more feasible in heavy-fermion super- 
conductors without inversion symmetry such as CePtsSi, 
where the enormous mass enhancement suppresses the 
magnitude of the Meissner supercurrent, than in conven- 
tional metals with moderate effective electron mass. We 
have proposed the experimental setup for the observation 
of the magnetoelectric effect in CePtaSi which utilizes the 
Volovik effect. It has been also pointed out that in non- 
centrosymmetric p wave superconductors, the coherence 
effect on the nuclear relaxation rate \JT\T is similar to 
that of conventional s wave superconductors. 

Finally, we would like to comment on the implication of 
our results for Ulr, which is the recently discovered ferro- 
magnetic superconductor without inversion symmetry.- 7 
Ulr exhibits superconductivity under high pressure in 



the vicinity of the phase boundary between ferromagnetic 
and non-magnetic states. The resistivity of this system 
increases remarkably as the applied pressure approaches 
to the critical value at which the ferromagnctism dis- 
appears, indicating the existence of ferromagnetic criti- 
cal fluctuation. In this case, the magnetoelectric coeffi- 
cient ea. (|36H may be enhanced by the three-point vertex 
functions A*?^ s, of which the magnitudes are much in- 
creased by ferromagnetic fluctuation, provided that the 
spin easy axis is taken as the x axis. Since, the crys- 
tal structure of Ulr is monoclinic, and does not possess 
any mirror planes, the Rashba spin-orbit interaction with 
the n vector perpendicular to the spin easy axis should 
always exists. Thus, the magnetoelectric effect strongly 
enhanced by ferromagnetic fluctuation may be observed 
in Ulr under an applied magnetic field parallel to the spin 
easy axis. 
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